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ABSTRACT 

Recent BOOMERANG and MAXIMA-I analyses of the angular power spectrum of the Cos- 
mic Microwave Background put more stringent constraints on the cosmological parameters. 
I show that these constraints are consistent with the observed baryon budget in clusters of 
galaxies, allowing any further contribution to this budget from baryons not in stars and not 
X-ray emitting to be less than 14 per cent at the 90 per cent level of confidence. 
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1 INTRODUCTION 

The calculations on the primordial nucleosynthesis abundance of 
the light elements (e.g. D, "^He, ''He, ^Li) give a direct measure- 
ment of the value of the baryon density, ilb, relative to the critical 
density necessary to close the Universe. If the regions that collapse 
to form rich clusters maintain the same ratio Jlb/f^c as the rest of 
the Universe, a measure of the cluster baryon fraction can be then 
used with the estimate of f2b to constrain the "cold", clustered com- 
ponent, Qc, of the total density parameter, Q,tot- 

The clusters baryon budget is composed mainly from the lu- 
minous baryons of the X-ray emitting gas that falls into the clus- 
ter dark matter halo. Other contributions come from the baryonic 
mass in galaxies, /gai, from intergalactic stars (Theuns & Warren 
1997) and a still hypothetical baryonic dark matter (like MACHOs, 
Alcock et al. 2000; but see Freese, Fields & Graff 2000). Given 
the large uncertainties on the relative contribution of the latter two 
sources to the baryon budget, I qualify these as "other baryons", 
fob - Therefore, one can put the following relation between the rel- 
ative amount of baryons in the Universe and in clusters with total 
gravitating mass, Mgrav, as inferred from the equation of the hy- 
drostatic equilibrium between the gas and the gravitational poten- 
tial that does not include the self-gravity of the hot plasma (e.g. 
Suto, Sasaki & Makino 1998): 



fib 



Mb 



/gas + /gal H~ fob > /g; 



(1) 



where fg, 



-Mgas /-A^grav , ./"gal — -^'^gal /^gTEiV , foh — 



Mob /Mg^^^ . (Note that any dependence upon the Hubble constant 
is discussed in the next section). In the following section, I discuss 
this equality in more details, considering the corrections required 
to compare the cluster baryon fraction to the universal value. 

When new tighter and lower constraints from nucleosynthesis 
were published by Walker et al. (1991; but see discussions in recent 
years on the abundance of the light elements, e.g. Hogan 1998), it 



became evident that just the amount of baryons visible in X-ray was 
enough to put in crisis an Einstein-de Sitter Universe with 0,^ = 1, 
giving rise to the so-called baryon catastrophe in clusters of galax- 
ies (White & Frenk 1991, White et al. 1993, White & Fabian 1995, 
David, Jones & Forman 1995, Evrard 1997, Ettori & Fabian 1999, 
Mohr et al. 1999; cf Fig. ^. 

Now that independent and well-constrained measurements of 
the cosmological parameters from estimates of the angular power 
spectrum of the Cosmic Microwave Background (CMB) are avail- 
able, one can reverse the problem and, using the estimated values 
for fib and f2c and the observed gas mass fraction, investigate the 
composition of the cluster baryon budget. 

In particular, one can estimate the ratio 



C = 



S^b/fic 



(2) 



that is expected to be close, but less than, the unity from eqn. |^. 

To do this, I use the estimates on the gas mass fraction pub- 
lished in Evrard 1997 (E97, as compilation of the measurements 
in White & Fabian 1995 and David, Jones, Forman 1995), Ettori & 
Fabian 1999 (EF99, for a sample of highly X-ray luminous clusters, 
iboi > 10*^ erg s"^), Mohr et al. 1999 (MME99, for an X-ray flux 
limited sample). I use the measurements quoted at rsoo, where the 
mean density in the cluster is 500 times the background value. The 
overdensity of 500 represents a confident outer limit where our as- 
sumptions on the intracluster gas as isothermal and in hydrostatic 
equilibrium still hold. 

Moreover, one can also include the stellar contribution to 
the baryon clusters budget. This is /gai — Mgai /Mgrav ~ 
0.020toml'^~^ (White et al. 1993, Fukugita et al. 1998). 

Hereafter, I refer to Qm as the total matter density in unity of 
the critical density, pc = 3Hq /(SivG), where Ho is the Hubble 
constant that we represent as Ho = 50/i5o km s~' Mpc~^. Hence, 
the total energy in the Universe in units of the critical value can be 
written as 
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2 



(fib + fic) + fiA + fik = fim + fiA + fik = 1 



(3) 



where S7a is the constant energy density associated with the "vac- 
uum" (e.g. Carroll et al. 1992) and fik is the energy density related 
to the curvature^. For a flat Universe, fik = and Sim + SIa = 1- 



2 THE COSMOLOGICAL CONSTRAINTS FROM 
BOOMERANG AND MAXIMA-I 

BOOMERANG (de Bemardis et al. 2000, Lange et al. 2000) and 
MAXIMA-I (Hanany et al. 2000, Balbi et al. 2000) are two baloon 
experiments that looking to the CMB on scales of few tens of ar- 
cminutes have permitted to constrain the cosmological models with 
a fit to the angular power spectrum of the temperature anisotropy in 
the detected signal. Jaffe et al. (2001) have combined these limits 
with the ones obtained from COBE/DMR (Bennett et al. 1996) on 
scales of few degrees. 

I consider here the best-fit results on Slb^iio = '^b and 
^ch^Q = i^c quoted in Table 1 of Lange et al. (2001) and the 
values for an assumed Sltot = 1 of ojb = 0.030 ± 0.004, u>c = 
0.19 ± 0.07, i/o = 75 ± 10 km s~^ Mpc"^ in Jaffe et al. (2001). 

The ratio cu^ l^c ~ flb/^c represents the relative amount of 
baryons with respect to the cold dark matter as measured in clusters 
of galaxies from the ratio of the baryonic mass to the total mass in 
hydrostatic equilibrium with the gas itself at a given radius. Using 
the gas fraction estimates in E97, EF99, MME99, I compute the 
ratio C that is expected to be < 1. 

Before to proceed to the measure of C, we have to correct the 
estimates of the observed /gas , that are generally expressed in terms 
of /i5o = 1 and fim ~ 1, for changes in the cosmological models 
and known systematic errors. These corrections are listed below 
and their contribution indicated as relative error, c,; = (/ga"'^ — 

/obs \ / robs 
gas // J gas ■ 

The gas mass fraction, /gas is given by the ratio of the gas 
mass and of the total gravitating mass within some radius, rsoo in 
the present study: 



Mgas(< J"500) 
Mgrav(< r-50o) 



{r)r^dr 



G^mp 



( 9 In pgaa \ 

V 81-- J, 



(4) 



where I assume that the plasma is isothermal and in hydrostatic 
equilibrium with the gravitational potential, r = 6 dang is the phys- 
ical radius and depends upon the angular separation, 0, and the an- 
gular diameter distance, dang, Tgas is the cluster gas temperature, 
yu is the mean molecular weight in a.m.u. (~ 0.6), G is the gravita- 
tional constant and rup is the proton mass. 

Equation ^ depends on cosmology to define (i) the physical 
radius, r, and (ii) the overdensity A here assumed equal to 500 for 
an "f2m = 1" Universe. 

The first cosmological correction, ci, is evaluated as follows. 
Given that the surface brightness, S{9), is the integration along the 
line of sight of the emissivity of the intracluster plasma due to ther- 
mal bremsstrahlung [i.e., S{d) ~ J p^^sTg^sdr], the gas density is 
proportional to cirng^- Combining this with the other dependence 
in eqn. W, 



(5) 



Contribution from neutrinos is not considered because their role, i.e. their 
mass, is still uncertain 



N clusters: 30 
Prob (all wrong)= 0.003% 
6B% range = [0.138, 0.204] 
90% range=[0.114, 0.232] 
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Figure 1. Bayesian probability distribution for the gas fraction observed in 
30 high-luminosity clusters (from Ettori & Fabian 1999). The distribution 
peaks at /gas = 0.168. The arrow indicates the baryon density of 0.076 
h"'^ (Buries & Tytler 1998) estimated through the primordial nucleosyn- 
thesis calculations from the observed low deuterium abundance. It has a 
probability of 7.2 X 10"'' with respect to the plotted distribution, making 
unlikely the agreement between the two baryon estimates in an Einstein-de 
Sitter Universe. 



I compute the angular diameter distance as a function of the 
cosmological parameters S7m + S7A+fik = 1 (cf. eqn. 25 in Carroll, 
Press & Turner 1992): 



da 



(6) 



where S{lu) is sinh((jj), lu, sin((jj) for fik greater than, equal to and 
less than 0, respectively, and 



1/2 



dC 



For f^k 



dans — 



[(l + C)2(l + fin.C)-C(2 + C)nA]'/'' 

0, eqn. ^ can be written as 

( z if flm — 



Ho{l + z) 



Jo 



-1/2 



[{i-i-C)3+n; 



11/2 



Otherwise 



(7) 



(8) 



The correction, ci, introduced by a low matter density Uni- 
verse and a Hubble constant h^o ~ 1-5 lower the gas fraction by 
about 40 per cent (ci ~ —0.40) mainly due to the dependence 
upon the Hubble constant (/gas oc /i^q^'^). As shown in Fig. ^ a 
low density Universe rises the gas fraction measurement for nearby 
clusters by about 10 per cent. 

The second cosmological correction, C2, appears as variation 
of the radius, ta, at which the given mean overdensity of the total 
mass within the cluster. A, with respect to the background value, 
fimPc is reached in different cosmological scenarios. To estimate 
A in function of the cosmological parameters, I use the relation 
valid for a "f2ni + ^a ~ 1" Universe (e.g. Kitayama & Suto 1996, 
Henry 2000; see Fig. |): 
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Figure 2. (Left panel) Ratio between the dependence of the gas fraction on the angular diameter distance varying f2ni and with f!m = 1- The same Hubble 
constant is assumed. (Right panel) Ratio between the overdensity, A, as function of redshift and Qm to the present value of 500 adopted in this study. 



A(an,2) = A(an = l) 



1 + 0.4093 (fi" 



-1 



,(9) 



where - 1) = (^T/ - 1)/(1 + zf. 

The measurement of ta requires an observed mass profile 
Mgrav(r), SO that Mgrav(< ta) = (4/3)7rriAanPc(l + 2:)^ 
Considering that Mgrav (r) oc r in the region of interest, 



rA 



■ r-500 



V 500 



(10) 



where 500 is the mean overdensity adopted in the samples of data 
considered here for an "fim = 1" Universe. Observations (David 
et al. 1995, EF99) and hydrodynamical simulations (see E97) agree 
on the radial dependence of the gas fraction near rsoo, /gas oc r"'^. 
Thus, the gas fraction has to be corrected by 



/gas(rA) = /gas(r500) (-^^) 



(11) 



This correction increases the gas fraction value by about 6 per cent 
or less (C2 ~ +0.06). 

A third correction, C3, is required from the evidence in both 
simulations and observations of baryon depletion in galaxy clus- 
ters. The combined N-body and gas dynamics simulations suggest 
that the amount of the cosmic baryons that falls in clusters as gas 
component at rsoo is slightly lower than 1 (~ 0.9, Frenk et al. 
1999; ~ 0.88 in MME99), in agreement with the observed increase 
of the gas fraction with radius as r"'^ near to rsoo (EF99, David 
et al. 1995). This implies that the estimated gas fraction (and the 
derived value C) might be underestimated by about 12 per cent 
(C3 « +0.12). 

A fourth correction, C4, comes still from hydrodynamical sim- 
ulations of galaxy clusters that indicate an overestimate of the gas 
mass due to the neglected clumping of the plasma. Considering 
that the observed thermal emission scales proportionally to Pgas> 
the presence of X-ray emitting clumps induce an overestimate on 
the gas fraction by about 12 per cent respect to the "true" value in 
simulated clusters, as shown by MME99 (C4 ~ —0.12). 



Other systematic errors present in the measurements of the 
gas fraction and whose magnitude is still uncertain are discussed in 
Section 2.2. 

Our analysis now requires an estimate of the "true" gas mass 
fraction given the four corrections listed above: 



/gas — /ga 



(12) 



Note that just the two ones related to the cosmological models as- 
sumed play a relevant role, because the third and fourth ones cancel 
(i.e. C3 + C4 ~ 0). 

I calculate a weighted mean value of the distribution in the 
gas fraction measurements (after the cosmological corrections are 
applied and including the quoted errors) for each dataset indepen- 
dently: 

r 0.108 + 0.005 (E97) 
/gas = < 0.115 + 0.002 (EF99) (13) 
[ 0.130 + 0.001 (MME99) 

I plot the results for C in Fig. ^ 

I observe that (i) C 1 for all the sets of cosmological val- 
ues discussed in Table 1 of Lange and collaborators (see Fig. |3|), 
providing evidence that there is now agreement between the fa- 
vorite cosmological scenario and the observed amount of baryons 
in galaxy clusters, (ii) none of these cosmological sets can be ex- 
cluded on the basis of the cluster baryon budget due to the large (a 
relative error of 40 per cent at la level) uncertainty still presents in 
the estimate of cun- 



2.1 Cluster baryon budget: gas, galaxies and ... then? 

From the estimate of C, and including the contribution from stars 
in galaxies, we are left with plausible constraints on any additional 
baryonic components, like baryons in form of dark matter or inter- 
galactic stars: 
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Figure 3. Left panel Estimates of the ratio C for gas fraction data in Evrard (1997, E97), Ettori & Fabian (1999, EF99), Mohr et al. (1999, MME99) and the 
joint cosmological constraints from BOOMERANG, MAXIMA-I and COBE/DMR in Jaffe et al. (2001). Right panel Assuming the gas fraction value from 
the compilation in Evrard (1997), I calculate C for the coiresponding hypotheses Pi (with i from to 1 1) in Lange et al. (2001) and for the best-fit values from 
the MAXIMA experiment. The diamonds indicate the central value given the only gas fraction, whereas the triangles include also the stellar contribution. The 
error bars come from the propagation of the uncertainties on the considered parameters at Icr level and refer to the estimates of C with the gas fraction only. 



/ob = - .fga. - ,f,al = (1 - C - X (14) 

From the best-fit results in Jaffe et al. (2001), I obtain that 
(error at Icr level): 

r 0.037 ± 0.064 (E97) 
/ob = < 0.030 ± 0.064 (EF99) (15) 
i 0.014 ± 0.064 (MME99) 

Therefore, no more than 14 per cent (90 per cent confidence level 
= 1.64(j; 23 per cent at 3cr = 99.7 per cent level of confidence) of 
the total matter in clusters of galaxies can be baryons not accounted 
for in X-ray emitting plasma and stars in galaxies. 



2.2 Systematic uncertainties on /gas 

The assumption of an isothermal gas in hydrostatic equilibrium 
with the gravitational potential provides an useful framework to 
work out the contribution of the X-ray emitting plasma to the clus- 
ter baryon budget. However, several aspects of the physics of the 
intracluster medium are still under investigations and can introduce 
some systematics errors on the measure of the gas mass fraction. 
I list here the most significant of these, considering their effect on 
the estimates of the values of C and /ob . 

• While the gas mass is well constrained from present X-ray ob- 
servations, the total mass profile still represents a significant uncer- 
tainty in the the baryon fraction estimates. Particularly, the shape 
of cluster temperature profiles is still defined with contradictory 
results (Markevitch 1998, Irwin et al. 1999, White 2000). I have 
assumed here an isothermal gas. Now, I investigate the variation of 
/gas once a polytropic profile for the intracluster gas (Tgas oc p^^^ ) 
is considered. I model the changes in the total gravitating mass us- 
ing a /3— model (Cavaliere & Fusco-Femiano 1976, Ettori 2000) 
with parameters — 0.3Mpc, 13 — 2/3 and a polytropic index, 7, 
equal to 1.24 (Markevitch 1998). 

• The presence of a magnetic field in the cluster plasma, like 
the recent reported measurements of Faraday Rotation (Clarke et 
al. 1999, Taylor et al. 1999) and excess in emission in the Hard 



X-ray band (Fusco-Femiano et al. 1999) seem suggest, might sup- 
port a non-thermal component that (i) contributes to the total pres- 
sure and (ii) mimics the thermal emission by a contribution of 
about 10-25 per cent (Fusco-Femiano et al. 1999) that has to be 
accounted for once the cluster luminosity is recovered from the 
X-ray count rate. Assuming a proportional dependence between 
the non- thermal pressure, Pnt, and the gas pressure, Pgas, such 
as Pnt = i^Pgas, the total pressure is then Ptot = (1 + !^)Pgas 
and affects the total gravitating mass through the hydrostatic equa- 
tion. I also include the correction by a factor e ~ \/0.80 ~ 0.90 
of the gas density due to the non-thermal contribution to the total 
observed emissivity that can be reduced . 

Referring to the corrected quantities as NEW with respect to 
the previous, OLD ones, I can write: 

Mgas.NEW ~ eAfgas.OLD 

Mgrav.NEW = (1 + 1^)7 A^grav, OLD (16) 

f g Tq r 

./gas, NEW — ^(1 + ^) Tgaa(r) .'giis.OLD 

I plot in Fig. ^ how the gas fraction estimated under the 
isothermal assumption changes for the corrections discussed here. 
I note that the presence of a temperature gradient reduces the to- 
tal mass and, therefore, increases the gas fraction value, whereas a 
non-thermal spectral component lowers the observed thermal con- 
tribution and, consequently, the estimate of the gas fraction at rsoo • 

Then, a negative gradient in the plasma temperature profile 
can reduce up to the 40 per cent the room available for any fur- 
ther baryonic contribution, only partially compensated (if below 
the limit of the equipartition with thermal plasma) by the presence 
of a non-thermal component. 



3 CONCLUSIONS 

Comparing the recent cosmological constraints from measurements 
of the angular power spectrum of the temperature anisotropy in the 
CMB with the observed distribution of the gas mass fraction in 
clusters of galaxies, I conclude that 
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Figure 4. Main uncertainties on the gas fraction represented as clianges 
with respect to /g!,° estimated from a plasma assumed isothermal and with 
no contribution from non-thermal component. The dashed line represents 
the gas fraction corrected by the presence of a temperature gradient with 
polytropic index, 7 = 1.24. The solid thick line includes a correction for a 
non-thermal component with the quoted e and v. The thin solid line assumes 
Ptot = 2Pgas. A magnetic field, B, of ~ fiG has an associated Pti'j: 
of about 50 K cm"'' ~ 10~'^Pgas. Note that the presence of a tempera- 
ture gradient tends to increase the gas fraction value at r ~ 5rc ~ f'soo, 
whereas correcting by the contamination of a non-thermal component re- 
duces the gas fraction. 



(i) the ratio C = /gas/(fib/fic) is close to 1 (Fig. |) for all 
the sets of cosmological values discussed in Jaffe et al. (2001) and 
Lange et al. (2001). Although I can not ruled out any of these mod- 
els just from cluster baryon budget due to the large relative error 
that affects ui^, on the other hand it shows that the most favorite 
cosmological scenario is consistent with the observed baryon bud- 
get in clusters of galaxies; 

(ii) there is evidence that only 14 per cent (90 per cent confi- 
dence level; 23 per cent at 3a level of confidence) of the total mass 
can be present in galaxy clusters in the form of baryons neither in 
stars nor X-ray emitting. It is worth to note that using the calcula- 
tions in Freese et al. (2000) and B-band luminosity and X-ray mass 
of the Coma Cluster in White et al. (1993), I estimate a contribution 
of MACHOs to the cluster mass, /macho = A^macho /Afgrav, of 
about 2 per cent, that is within the above constraint. 

(iii) the upper limit of 14 per cent is in the order of the observed 
scatter in the gas mass distribution of about 15-20 per cent (cf. 
Fig. |l|), suggesting that part of this scatter might be due to still un- 
observed baryons. It is worth to note, however, that the present ac- 
curacy of the gas mass measurements is in the order of 10 per cent 
(cf. simulations in MME99) and probably larger uncertainty is still 
related to the gravitating mass estimates (see next item). The com- 
bination of these uncertainties could explain the observed scatter in 
the gas mass fraction distribution. Data from Chandra and XMM- 
Newton X-ray observatories are in condition to reduce significantly 
the size of these uncertainties; 



(iv) systematic errors on the gas mass fraction measurement, 
like the presence of a temperature gradient in the plasma or a non- 
thermal component in the cluster emission, tends to reduce by a 
value between 25 and 40 per cent the estimate of the fraction of the 
gravitating mass as "unseen" baryons. 
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